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Abstract 

In this work we further extend the investigation of holographic gauge theories in external 
magnetic fields, continuing earlier work. Wc study the phenomenon of magnetic catalysis 
of mass generation in 1+3 and 1+2 dimensions, using D3/D7- and D3/D5-brane systems, 
respectively. We obtain the low energy effective actions of the corresponding pseudo 
Goldstone bosons and study their dispersion relations. The D3/D7 system exhibits the 
usual Gell-Mann-Oakes-Renner (GMOR) relation and a relativistic dispersion relation, 
while the D3/D5 system exhibits a quadratic non-relativistic dispersion relation and a 
modified linear GMOR relation. The low energy effective action of the D3/D5 system is 
related to that describing magnon excitations in a ferromagnet. We also study properties 
of general Dp/Dq systems in an external magnetic field and verify the universality of the 
magnetic catalysis of dynamical symmetry breaking. 
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1 Introduction 



The applications of holographic gauge/gravity correspondences to the study of more and more 
diverse phenomena are ever widening in scope. Over the last half a decade the links between 
finite temperature generalizations of AdS/CFT and experimental heavy-ion collisions have be- 
come much more concrete and the theoretical methods available to us are yielding ever deeper 
results concerning the properties of the quark-gluon plasma (see ref. [HIS] for a recent review). 
Moreover in the last year the links between holography and condensed matter systems have 
also flourished, with work on superconductivity, and superfluidity, quantum phase transitions 
and both the classical and quantum hall effects having recent successes (e.g. refs. [31 H] and 
references therein). 

In the present work we extend the investigation of holographic gauge theories in the presence of 
external magnetic flelds from the work flrst studied in ref. [5]. In this paper we are interested in 
finding both universal properties of strongly coupled gauge theories in the presence of magnetic 
fields, as well as in the different phenomena exhibited in such theories in a variety of space-time 
dimensions. 

The phenomenon of dynamical flavor symmetry breaking catalysed by an arbitrarily weak 
magnetic fleld is known from refs. [6l [7] and refs. [8l |9l |10]. This effect was shown to be model 
independent and therefore insensitive to the microscopic physics underlying the low energy 
effective theory. In particular the infra-red (IR) description of the Goldstone modes associated 
with the dynamically broken symmetry should be universal. We therefore expect to be able 
to study this phenomenon using the holographic formalism. The aim of the present study will 
be to investigate the dynamics of the Goldstone modes and construct the low energy chiral 
Lagrangian of theories both in 3+1 and 2+1 dimensions in the presence of external magnetic 
fields, showing that the appropriate holographic models give precisely the results expected from 
the traditional field theory approach. 

The effective dynamics of fermion pairing, in c? + 1 dimensions, in the presence of an external 
magnetic field is constrained to d — 2 spatial dimensions. For this reason there are marked 
differences in the phenomenology of such systems in two and three spatial dimensions. In 2+1 
dimensions refs. [7j-[in] Poincare symmetry is broken by the magnetic field (there is no longer 
any trace of the original boost invariance), removing the strong constraints on the dynam- 
ics of Goldstone modes imposed by special relativity. The naive Goldstone boson counting 
therefore does not hold and the resulting dispersion relation for the Goldstone modes takes a 
quadratic form, unlike in the case of 3+1 dimensions where an S0{1, 1) subgroup of 5*0(3, 1) 
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constrains the dynamics. Although the number of Goldstone particles is no longer constrained 
in the non-relativistic setting, the number of Goldstone fields is fixed by the dimension of G/H 
((j>'=symmetry of the action, if=symmetry of the ground state). We will show that this also 
holds in the AdS/CFT context. 

The 2+1 dimensional model is of particular interest because, as shown in ref. [TT], the low 
energy effective description is that of magnon excitations in a ferromagnet. Using a D3/D5 
brane intersection we will be able to reproduce such a result at quadratic order in the chiral 
Lagrangian. Moreover such 2+1 dimensional theories may have relevance in the arenas of the 
quantum hall effect, and high Tc superconductivity. 

In addition to the phenomena discussed specifically in two and three spatial dimensions we show 
that certain universal behaviors are exhibited holographically in the present context. Here we 
will study holographic systems T-dual to the D3/D7 flavor model and show that the existence 
of an arbitrarily small magnetic field induces a spiral behaviour in the equation of state for 
such systems. In the limit that the chiral symmetry of the underlying theory is preserved, this 
equation of state can be studied analytically and such a symmetric vacuum can be shown to be 
unstable. This is holographically equivalent to the findings of refs. [6]- [10] - fiavor symmetry 
breaking is induced dynamically by the presence of a magnetic field. 

The outline of the present paper is as follows: 

In section 2 we will return to the D3/D7 brane intersection in the presence of an external 
magnetic field, discussed in ref. [5j. We shall show explicitly how the magnetic catalysis of 
flavor symmetry breaking is realised in the holographic system, including the calculation of the 
chiral Lagrangian to second order in the low energy degrees of freedom. We will show that the 
Gell-Mann-Oakes-Renner relation holds analytically and obtain the dispersion relation for the 
Goldstone modes. 

In section 3 we turn to the case of the D3/D5 defect theory and show here how the 5*0(3) 
flavor symmetry is dynamically broken to the U{1) subgroup in the presence of a magnetic 
field. In this non-relativistic system we find the Goldstone modes and show that the number of 
massless modes is not the same as the number of broken generators, but satisfies a more general 
counting rule [25] applicable for non-relativistic systems. We also show that the single Goldstone 
mode satisfies a modified Gell-Mann-Oakes-Renner relation and a quadratic dispersion relation. 
Again we obtain the dispersion relation analytically in the small mass limit and find the low 
energy effective Lagrangian which describes magnon excitations in a ferromagnet. 

In section 4 we prove that the magnetic catalysis of dynamical symmetry breaking is a universal 
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effect in gauge theories dual to Dp/Dq intersections in the appropriate decouphng hmit. This 
proof involves showing that all such systems exhibit a self-similar spiral behaviour in their 
equation of state which leads to an instability for the solution with zero dynamical mass. Just 
as in the work of ref. [7j this effect is independent of the magnitude of the external magnetic 
field. 

2 Mass generation in the D3/D7 system 

In this section we will review the results of refs. [HI [121 [13], where a holographic study of flavored 
A/" = 4 supersymmetric Yang- Mills in an external magnetic field was studied using the D3/D7 
system. We will focus on the effect of mass generation by magnetic catalysis in this theory 
and provide a detailed analysis of the pseudo-Goldstone mode associated to the spontaneous 
breaking of a global U{1) R-symmetry. In particular we will show that the Gell-Mann-Oakes- 
Renner relation for the mass of the corresponding t]' meson is satisfied. 

The D3/D7 system provides a dual holographic description of Nf fundamental M = 2 hy- 
permultiplets coupled to A/" = 4 SU{Nc) supersymmetric Yang Mills theory in the quenched 
approximation Nf <^ [14J. At zero separation between the D3 and D7-branes the funda- 
mental hypermultiplets are massless and the /^-function of the theory is proportional to Nf/Nc. 
Thus in the quenched approximation the /3-function vanishes and the corresponding gauge the- 
ory is conformal. The global 5*0(6) R-symmetry of the A/" = 4 SYM theory is broken to an 
SU(2) X f/(l) R-symmetry, the U{1) corresponding to rotations in the 2-plane transverse to 
both the D3 and D7-branes. The left and right handed fermions of the hypermultiplet have 
opposite charges under this U{1)r and thus the formation of a fermionic condensate (ipip) would 
lead to the spontaneous breaking of this symmetry. 

2.1 Spontaneous symmetry breaking 

There are various ways in which one can study the breaking of the chiral symmetry holo- 
graphically. This has been studied in the past by the deformation of AdSs x S*^ by a field 
corresponding to a marginally irrelevant operator on the gauge theory side refs. [T71 HBl [12]. 
In the present case however we will stimulate the formation of a condensate by turning on the 
magnetic components of the f/(l) gauge field of the D7-branes F^/j (equivalent to exciting a 
pure gauge S— field in the supergravity background). This U{1) gauge field corresponds to 
the diagonal U{1) of the full U{Nf) gauge symmetry of the stack of D7-branes. Since the 
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D7-branes wrap an infinite internal volume, the dynamics of the U{Nf) gauge field is frozen in 
the four dimensional theory and the U{Nf) gauge symmetry becomes a global flavor symmetry 
U{Nf) = U{1)b X SU{Nf). Therefore the U{1) gauge field that we consider corresponds to 
the gauged U{1)b baryon symmetry and the magnetic field that we introduce couples to the 
baryon charge of the fundamental fields |27j . 

The problem thus boils down to studying embeddings of probe D7-branes in the AdSs x 
background parameterized as follows: 

ds^ = ^-—^ — [-dxl + dxl + dxl + dxl] + ^ [dp^ + p^d^l + dL^ + L^dcj)^] , 

dnj = di:^ + cos^ i:dp^ + sin^ ijd'j^ , (1) 



where p, if), jS, 7 and L, (j) are polar coordinates in the transverse and planes respectively. 

Here Xa=i..3, p, i', P, 7 parameterize the world volume of the D7-brane and the following ansatz 
is considered for its embedding: 

= const , L = L{p) , 
leading to the following induced metric on its worldvolume: 

d~s = P-^^^ [-dxl + dxl + dxl + dx'^ + -^r^^, [(1 + L\pf)dp' + p'd^ll] . (2) 

The D7-brane probe is described by the DBI action: 

^DBi = -iV//i7 j c/'ee-*[-det(G„, + Bab + 2Tia' Fab)]^'^ . (3) 
Ma 

Here pi = [(27r) ''a;'"^]^^ is the D7-brane tension, Gab and Bab are the induced metric and i?-field 
on the D7-brane's world volume, while Fab is its world-volume gauge field. A simple way to 
introduce a magnetic field is to consider a pure gauge i?-field along the X2, X3 directions: 

S(^) = Hdx2 A dx3 . (4) 

Since Bab and Fab appear on equal footing in the DBI action, the introduction of such a S-field 
is equivalent to introducing an external magnetic field of magnitude H/ {2na') to the dual gauge 
theory. 
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Though the full solution of the embedding can only be calculated numerically, the large p 
behaviour (equivalently the ultraviolet (UV) regime in the gauge theory language) can be 
extracted analytically: 



L{p) 



m 



+ 



(5) 



As discussed in ref. [TH], the parameters m (the asymptotic separation of the D7- and D3- 
branes) and c (the degree of bending of the D7-brane in the large p region) are related to 
the bare quark mass rriq = m/2'Ka' and the fermionic condensate {ipip) oc — c respectively. It 
should be noted that the boundary behavior of L(r) really plays the role of source and vacuum 
expectation value (vev) for the full M = 2 hypermultiplet of operators. In the present case, 
where supersymmetry is broken by the gauge field configuration, we are only interested in the 
fermionic bilinears and this will refer only to quarks, and not their supersymmetric counterparts. 

At this point it is convenient to introduce dimensionless parameters c = c/B?H^/^ and m = 
m/R\fH. By performing a numerical shooting method from the infrared while varying the 
small p boundary value, L{p 0) = Lir, we recover the parametric plot presented in figure [l| 
the main result explored in ref. [H]. 
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Figure 1: Parametric plot of c against fa for fundamental matter in the presence of an external 
magnetic field. The lower (black) line represents the curve 1/m, fitting the large rh behavior. It is 
also evident that for the outer branch of the spiral, for m = the condensate, (V'V') is non-zero. The 
corresponding value of the condensate is Ccr = 0.226. 



The lower (black) curve corresponds to the analytic behavior of c(m) = 1/m for large fh. The 
most important observation is that at m = there is a non-zero fermionic condensate: 



H 



3/2 



(6) 
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Where A = Qym^c is the 't Hooft couphng and ~ 0.226 is a numerical constant corresponding 
to the y-intercept of the outer spiral from figure [l] Equation ^ is telling us that the theory 
has developed a negative condensate that scales as (2^)^^^- This is not surprising, since the 
theory is conformal in the absence of the scale introduced by the external magnetic field. The 

1 /2 

energy scale controlled by the magnetic field, {1^^) , leads to an energy density proportional 
to (^^)^- In order to lower the energy, the theory responds to the magnetic field by developing 
a negative fermionic condensate. 

Another interesting feature of the theory is the discrete-self-similar structure of the equation 
of state (c vs. m) in the vicinity of the trivial m = embedding, namely the origin of the plot 
from figure [T] presented in figure [2j 
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Figure 2: A magnification of figure [T] shows the spiral behavior near the origin of the (— c, ?fi)-plane. 
The second (left) spiral arm represents the (m, — c) — > (— m, c) symmetry of the theory. 

This double logarithmic structure has been analyzed in ref. [12], where a study of the meson 
spectrum revealed that only the outer branch of the spiral is tachyon free and corresponds to 
a stable phase having spontaneously broken chiral symmetry. In Section 3 of this paper we 
will show that an identical structure is also present for the D3/D5 system and in Section 4 
we will demonstrate that this structure is a universal feature of the magnetic catalysis of mass 
generation for gauge theories holographically dual to Dp/Dq intersections. 

A further result of refs. [5l [121 [13] was the detailed analysis of the light meson spectrum of 
the theory. In ref. [5] it was shown that the introduction of an external magnetic field breaks 
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the degeneracy of the spectrum studied in ref. [15]. This manifests itself as Zeeman sphtting 
of the energy levels. In the limit of zero quark mass, the study also revealed the existence 
of a massless meson" corresponding to the spontaneously broken f/(l)R symmetry. In the 
next subsection we will revisit the study of the meson spectrum of the theory and provide an 
analytic proof of the Gell-Mann-Oakes-Renner relation |24j : 

Ml = -^ji^m, , (7) 

in the spirit of the analysis performed in ref. |18j . 

2.2 The Gell-Mann— Oakes— Renner relation - an analytic derivation 

In order to study the light meson spectrum of the theory one needs to consider the quadratic 
fluctuations of the D7-brane embedding and study the corresponding normal modes [12]. Tech- 
nically one should consider the full supergravity action for the D7-branes: 

'S'tot = 5'dbi + 'S'wz , (8) 
where ^dbi is given by equation ([s]) and the relevant part of the Wess-Zumino term is given by 

Swz = ^^^^f^7 j F^2) A F(2) A C(4) + (27ra')/W7 j i"(2) A 5(2) A P[C(4)] , (9) 
The next step is to consider fluctuations of the D7-brane in the transverse M^: 

L = Lo{p) + 27ia'SL ; = 2na'<^ , (10) 

and expand equation ([s]) to second order in a'. Note that with such an expansion we should also 
consider fluctuations of the U{1) gauge field on the D7-brane. As demonstrated in refs. [5] [20] 
the effect of the magnetic field will be to mix the equations of motion for the scalar and vector 
fluctuations. In particular $ couples to the Aq and Ai components of the gauge field, while 6L 
couples to the A2 and ^3 components. The rest of the components of the vector field decouple 
and can be consistently set to zero. This splitting of the meson spectrum is a manifestation of 
the broken Lorentz symmetry. Indeed the external magnetic field breaks the 5*0(1,3) Lorentz 
symmetry down to 5*0(1, 1) x 50(2) corresponding to boosts in the Xq, Xi plane and rotations in 
the X2, X3 plane. Since the massless "pion" that we are interested in corresponds to fluctuations 
along 0, we will excite only the $, Aq, Ai fields. The relevant terms of the expansion are [5]: 
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= -{27rar^l^/M^^,S'^'^M^ , (11) 
Ua = -{27ra'y^y^\Hd,K^Foi , 



where: 



Gil Gil 

\S \\ = dia.g{—Gii ,Gii , +if2'(^2 _|_ ^pp ' ^i/'^/" ^"a' ^z?/?} ' (■^^) 



Here Lo(p) corresponds to the classical embedding of the D7-brane and Gab are the components 
of the background metric equation ([T]). 

The equations of motions for $ and Fqi are calculated from the quadratic action, resulting in: 



g{p) '\ 1 + ; p2 (^2 + ^(^) 
1 ^ /^(p)a,Foi\ , , i?^ Hd,K^ ^2 , n2^ 



9ip) "V 1 + ^? ; P' {P' + LIY 9ip) 
where Fqi = doAi — 9iAo and the gauge constraint —doAo + diAi = is imposed (note that 
this is the usual Lorentz gauge, corresponding to the unbroken S0{1, 1)) and we have defined: 

a = -d^, + d! + -^^^. (14) 



Once again the broken Lorentz symmetry is manifest in equation (14). The definition of the 
spectrum is now a subtle issue in the presence of the broken space-time symmetry. We will 
define the spectrum as the energy of a particle as measured in its rest frame. In fact because 
we retain the 5*0(1, 1) symmetry we may consider fluctuations propagating in the Xi direction. 
Since we are interested in describing the lowest lying modes ("pions" in particular) we will 
focus on modes that have no 5*^ dependence. Therefore we consider the ansatze: 

$ = e^(^°^"+^^^')/;,(p); Fqi = e*('=o^°+'=i^')/(p) , (15) 
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and define: 



M^ = kl-kl . (16) 



The equations ( 13 ) simplify to: 



g '\l + L'^ J (p2 + L^)2 g 

Note tfiat for large bare masses m (and correspondingly large values of L) the term proportional 
to the magnetic field is suppressed and the meson spectrum should approximate to the result 
for the pure AdSs x space-time case studied in ref. [15], where the authors obtained the 
following relation: 

2?77 / 

M„ = — v/(n + l)(n + 3) , (18) 

between the eigenvalue of the n^^ excited state uj^ and the bare mass m. If one imposes the 
boundary conditions: 

h{e) = l- /i'(e) = 0; /(e) = 1; /(e) = , (19) 

the coupled system of differential equations can be solved numerically. Then by requiring the 
functions h{p) and /(p) to be regular at infinity one can quantize the spectrum of the fluctu- 
ations. It is also convenient to define the following dimensionless parameter M = MR/\/H. 
The resulting plot for the first three excited states is presented in figure [3j There is Zeeman 
splitting of the states due to the magnetic field. (In the absence of the field there are three 
straight lines emanating from the origin; these are split to form six curves.) Also, at zero bare 
quark mass there is indeed a massless Goldstone mode, appearing at the end of the lowest 
curve. Furthermore the plot in figure |4] shows that for small bare quark mass one can observe 
a characteristic M oc \ffti dependence. In the next section we shall provide an analytic proof 
of that relation and obtain an integral expression for the numerical coefficient 0.64 presented 
above the plot in figure |4} 

In the following section we shall demonstrate that for small bare quark mass, irtq = m/2TTa', 
the spectrum exhibits the characteristic oc m dependence. Once we have illustrated that 
the functional dependence is correct we will show that the constant of proportionality is also 
that expected from the GMOR relation. Furthermore we shall generalize the ansatze (15) to 
consider fluctuations depending on both the momentum along the magnetic field k\\ = {ki, 0, 0) 
and the transverse momentum k± = (0, k2, k^): 

$ _ ^i{M.x)^^p^ . ^ ^i{ut+k.x)ji^p^ _ ^20) 
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Figure 3: There is Zeeman splitting of the states due to the magnetic field. In the absence of the 
field there are three straight lines emanating from the origin; these are split to form six curves. At 
zero bare quark mass (the end of the lowest curve) there is indeed a massless Goldstone mode. The 
straight lines correspond to the asymptotic AdS results. 

We shall also show that for small u = and \k\ the following dispersion relation holds: 

uj{kf = + kf\+-fkl ; uj = ko- A^ii = (A;i,0,0) ; ^± = (0,^2,^3), (21) 
where 7 is a constant that we shall determine. 



2.2.1 The oc m dependence 



Using an approach similar to the one employed in ref. [18] we define: 

gLl p4 1+^0 ~ - l + ^n 1 

U = It 



^0 \y ^ ^0) 



(22) 



The equations of motions (13) can then be written in the compact form: 



Let us remind the reader that for large p, Lq{p) has the behavior: 



(23) 



Lo oc m + — + 

p2 



(24) 
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Let us denote by Lq the classical embedding corresponding to (m = 0, c = Ccr). It is relatively- 
easy to verify that ai m = 0, kj_ = and correspondingly = uj"^ — kf, = the choice: 



ilj = ^ = ; ^1 = 



(25) 



is a solution to the system (23). Next we consider embeddings corresponding to a small bare 
quark mass 6m. This will correspond to small nonzero values of and k'j_. It is then natural 
to consider the following variations: 



^ = ^ + , 



(26) 



where 6ip and Sipi are of order M^. Note that M corresponds to the mass of the ground state 
at rriq = 6m/2TTa' and we are assuming that the variations of the wave functions Sip and 6ipi 



are infinitesimal for infinitesimal rriq. After expanding in equation (23) we get the following 
equations of motion: 



(27) 



where z/ = The second equation in (27) can be integrated to give: 

^i6ipi - ^i6i/ji = HK{uj^ - kl) + constant . 



(28) 
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From the boundary conditions that K\p=Q = and \l/i(0) = 0,\&i(0) = we see that the 
constant of integration is zero and arrive at: 



^1 



^2 



Next we multiply the first equation in (27) by \E' and integrate along p to obtain: 



(29) 



{u' - q) I dpj?^' - k'i I dpi>^'' = - I {^6^^J - ^6tlj)dp + I ^^6{^]dp+ (30) 



+H J ^^^^dp = -(^# - ^6ij) J + {^6^ - M^) 





5Jh 
^1 



dp , 



where the last term on the right-hand side of equation (30) has been integrated by parts 



using the fact that Sipi should be regular at infinity. From the definition of ^ it follows that 
^ oc p^/^Lo(0) as p — » and 'I' oc c/ p^^"^ as p ^ oo. This together with the requirement that 
ijji is regular at p = and vanishes at infinity, suggests that the first term on the right-hand 



side of equation (30) vanishes. For the next term, we use the fact that: 

1/4 ^ ^.^4 X 1/4 



1 + ^0^ 



12 



Lo + P' 



3/2 



12 



1 + ^'o 



and therefore obtain: 



5*|o = 0; 5^|o = 0, 

3 

(5^|oo oc p^l'^bm ; 64l\oo oc -y/p5m . 



The second term in equation (30) then becomes: 



2c6m . 



Finally using the equality in equation (29) we arrive at the result: 

oo _ oo 

2 l2\ / J„ J ,-,,T,2 , ^ ^ ( u2 j„r,,T,2 



Now we define: 



cu' - kf^) I dp { v^' + } -kl I dpv^' = 2c5m . 



oo \ / oo 

7 = I / dp^^^ \ / Ij dp{u^^ 



(31) 



(32) 



(33) 



(34) 



(35) 
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and solve for from equation (21) to obtain: 







2c6m . 



(36) 



Equation (36) suggests that the mass of the "pion" associated to the softly broken global U{1) 



symmetry satisfies the Gell-Mann-Oakes-Renner relation 

2(#) 



Mi 



2 



(37) 



In order to prove equation (37) we need to evaluate the effective coupling of the "pion" 
Noting that 5m oc and c oc —{ipip), we conclude that: 



(38) 



At this point is useful to verify the consistency of our analysis by comparing the coefficient in 



equation (36) to the numerically determined coefficient 0.64 from the plot in figure 4 Indeed 



from equation (36) we obtain: 

M/V^ -- 




-1/2 



0.655 



(39) 



where we have defined the dimensionless quantities: 

z> = H'^u; ^2 = m^/R^H^'^] ill = ml/R^H^'^] K = K/R"^ 
There is excellent agreement with the fit from figure |4j 



(40) 



Next we will obtain an effective four dimensional action for the "pion" and from this derive an 
exact expression for 

2.2.2 Effective chiral action and 

In this section we will reduce the eight dimensional world-volume action for the quadratic 
fiuctuations of the D7-brane to an effective action for the massless "pion" associated to the 
spontaneously broken global U{1) symmetry. Note that our effective action should be describing 



a single "pion" mode, while the 8D action given by equation (11) describes the dynamics of 



two independent degrees of freedom, namely $ and Fqi coupled by the magnetic B-field via 
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the second equation in equation (11 ). As rigid rotations along correspond to chiral rotations, 
(the asymptotic value of (f) at infinity corresponds to the phase of the condensate in the dual 
gauge theory) the spectrum of $ at zero quark mass contains the Goldstone mode that we are 
interested in. This is why we first integrate out the gauge field components Aq and Ai and 
then dimensionally reduce to four dimensions. 

Furthermore as mentioned earlier, because of the magnetic field the S0{1, 3) Lorentz symmetry 
is broken down to 5*0(1, 1) x 5*0(2) symmetry. This is why in order to extract the value of 
we consider excitations of $ depending only on the xq, xi directions and read off the coefficient 
in front of the kinetic term. The resulting on-shell effective action for $ is: 

S'''' = j d'x [-{do^f + (di^Y] , (41) 

where A/" is given by: 

oo _ 

(42) 



AT = (2W)2^Ar^7r2 y dp + ^ 



We refer the reader to Appendix A for a detailed derivation of the 4D effective action 5"^^. 

We have defined $ via cf) = (27ra;')$, where (f) corresponds to rotations in the transverse 
plane and is the angle of chiral rotation in the dual gauge theory. The chiral Lagrangian is 
then given by: 

= _(27ra')'^ J d^xd^^d''^^ ; /x = or 1 , (43) 



and therefore: 



oo 



The D7-brane charge in equation (44) is given by fi^ = [(27r)^a''*] ^ and the overall prefactor 
in equation (44) can be written as NfNc/2{2TTa')'^X. Now, recalling the expressions for the 



fermionic condensate, equation (|6]), and the bare quark mass, rUg = m/27ia', one can easily 



verify that equation (36) is indeed the Gell-Mann-Oakes-Renner relation: 



Ml = --^m, . (45) 



It turns out that for small momenta and small mass one can obtain the following 

more general effective 4D action (see appendix A for a detailed derivation): 

SeE = -U I rf^x|[-(9ol>)^ + (9il>)2]+7[(a2l>)^ + (93<l)V^^"^.^''}+--- , (46) 
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where 7 is defined in equation (132). As one can see, the action (46) is the most general 
quadratic action consistent with the 5*0(1, 1) x 5*0(2) symmetry and suggests that pseudo 
Goldstone bosons satisfy the dispersion relation ([2l]). 



3 Mass generation in the D3/D5 system 

In this section we provide a holographic description of the magnetic catalysis of chiral symmetry 
breaking in 1 + 3 dimensional SU{Nc) A/" = 4 supersymmetric Yang-Mills theory coupled 
to Nf M = 2 fundamental hypermultiplets confined to a 1 + 2 dimensional defect. Recently 
this theory received a great deal of attention and emphasis has been made of the potential 
application of this brane configuration in describing qualitative properties of 1 + 2 dimensional 
condensed matter systems (see for example refs. [281 EHl |30]). In this section we will study the 
effect of an external magnetic field on the theory and demonstrate that the system develops 
a dynamically generated mass and negative fermionic condensate leading to a spontaneous 
breaking of a global 50(3) symmetry down to a f/(l) symmetry. On the gravity side this 
symmetry corresponds to the rotational symmetry in the transverse M? . Naively there should 
be two massless Goldstone bosons corresponding to the generators of the coset 5'0(3)/f/(l). 
As we will show the 1+2 dimensional nature of the defect theory leads to a coupling of the 
transverse scalars corresponding to the coset generators and as a result there is only a single 
Goldstone mode. Furthermore the characteristic M^^ oc ^/m Gell-Mann-Oakes-Renner relation 
is modified to a linear M^^ oc m behavior. It turns out that these features can be understood 
from a low energy effective theory point of view. Indeed in 1 + 2 dimensions the effect of the 
magnetic field is to break the 5*0(1,2) Lorentz symmetry down to 5*0(2) rotational symmetry 
and as a result the theory is non-relativistic. A single time derivative chemical potential term 
is allowed (there is no boost symmetry) and interestingly the supergravity action generates 
such a term through the Wess-Zumino contribution of the D5-brane. It is this term that is 
responsible for the modified counting rule of the number of Goldstone bosons |25] and leads 
to a quadratic dispersion relation as well as to the modified linear Gell-Mann-Oakes-Renner 
relation. Another interesting feature of the model is that to quadratic order the effective low 
energy action is the same as the effective action describing spin waves in a ferromagnet [H] in 
an external magnetic field. We comment briefly on the possible applications of this similarity. 
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3.1 Generalities 



Let us consider the AdSs x supergravity background ([T]) and introduce the following param- 
eterization: 

ds^ = [-dxl + dxl + dxl + dxg] + [dr^ + r'^d^l + df + fd^'^^ , (47) 
= + ; dUl = da^ + cos^ ; dClj = + cos^ V'd^^ . 

We have split the transverse to x and introduced spherical coordinates r, and /, ^2 
in the first and second planes respectively. Next we introduce a stack of probe Nf D5-branes 
extended along the xo,xi,X2 directions, and filling the part of the geometry parameterized 
by r,fl2- As mentioned above on the gauge theory side this corresponds to introducing Nf 
fundamental M = 2 hypermultiplets confined on a 1 + 2 dimensional defect. The asymptotic 
separation of the D3 and D5 -branes in the transverse space parameterized by I corresponds 
to the mass of the hypermultiplet. In the following we will consider the following anzatz for a 
single D5-brane: 

/ = /(r) ; = ; = 0. (48) 

The asymptotic separation m = /(oo) is related to the bare mass of the fundamental fields via 
iriq = m/2Tra'. If the D3 and D5 branes overlap, the fundamental fields in the gauge theory are 
massless and the theory has a global S0{3) x S0{3) symmetry. Clearly a non-trivial profile of 
the D5-brane l{r) in the transverse would break the global symmetry down to S0{3) x f/(l), 
where U{1) is the little group in the transverse M^. If the asymptotic position of the D5-brane 
vanishes (m = 0) this would correspond to a spontaneous symmetry breaking, the non-zero 
separation /(O) on the other hand would naturally be interpreted as the dynamically generated 
mass of the theory. 

Note that the D3 /D5 intersection is T-dual to the D3 /D7 intersection from the previous section 
and thus the system is supersymmetric. The D3 and D5 -branes are BPS objects and there is 
no attractive potential for the D5-brane, hence the D5-brane has a trivial profile / = const. 
However a non-zero magnetic field will break the supersymmetry and as we are going to demon- 
strate, the D5-brane will feel an effective repulsive potential that will lead to dynamical mass 
generation. In order to introduce a magnetic field perpendicular to the plane of the defect, we 
consider a pure gauge 5-field in the xi, X2 plane given by: 

B = Hdxi A dx2 . (49) 

This is equivalent to turning on a non-zero value for the 0, 1 component of the gauge field on 
the D5-brane. The magnetic field introduced into the dual gauge theory in this way has a 
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magnitude H/2Txa' . The D5-brane embedding is determined by the DBI action: 

^DBi = -iV/yUs j rf'ee"*[-det(G,6 + Bab + 2iia'Fab)]"^ . (50) 

Me 

Where Gab and Bab are the pull-back of the metric and the S-field respectively and Fab is the 
gauge field on the D5-brane. 



With the anzatz (48) the Lagrangian is given by: 



/:ocr^v^rT7^^1 + p^. (51) 
From this it is trivial to solve the equation of motion for /(r) numerically, imposing Z(0) = Un 



and /'(O) as initial conditions. Clearly, at large r the Lagrangian (51) asymptotes to that at 
zero magnetic field and hence we have the asymptotic solution |2Bj: 

/(r) = m + ^ H , (52) 

where c oc ('?/''?/') the condensate of the fundamental fields. 

3.2 Spontaneous symmetry breaking 

Before solving the equation of motion it is convenient to introduce dimensionless variables: 

f = r/Ry/H ; I = 1/rVh ; m = tti/RVh ; d = c/R^H . (53) 
The Lagrangian ([5T|) can then be written as: 



C(xr^\/l + l'^Jl + . (54) 



The corresponding equation of motion is: 



Before solving equation (55) it will be useful to extract the asymptotic behavior of c(m) at 
large m. To this end we use that at large m the separation /(f) ^ m = const. The equation of 
motion then simplifies to: 

(r^ + m'^j-^ 
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and hence: 



-2m / dr 



(r^ _|_ 777,2)2 



(57) 



Using the expansion (52) one can verify that: 



hm r I' 

•— >+oo 



2m / dr 







im"^ 



(58) 



Equation (58) can thus be used as a check of the accuracy of our numerical results. Indeed the 

-c vs. fh is presented in figure |5} The most important observation 



numerically generated plot of 
is that at zero bare mass fh the theory has developed a negative condensate {ipip) oc — Ccr ~ 
—0.59. It can also be seen that for large fh the numerically generated plot is in good agreement 



with equation (58) represented by the lower (black) curve. Another interesting feature of the 



equation of state is the spiral structure near the origin of the parameter space analogous to 
the one presented in figure |2] for the case of the D3/D7 system. We will come back to this in 
Section 4 in more general terms, and show that this feature is universal for the class of gauge 
theories dual to the Dp/Dq systems. 




Figure 5: A plot of — c vs. rh. At zero bare mass m = the theory has developed a negative 
condensate (V'V') 

oc — Ccr ~ — 0.59. For large tti there is excellent agreement with equation 058^, as 
represented by the lower (blue) curve. 

In order to show that the global 5*0(3) symmetry is indeed spontaneously broken we need to 
study the free energy of the theory. Indeed the existence of the spiral structure suggests that 
there is more than one phase at zero bare mass, corresponding to the different y-intercepts of 
the — c vs. fh plot. We will demonstrate below that the lowest positive branch of the curve 
presented in figure |5] is the stable one. 
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Following ref. [26] we will identify the regularized wick rotated on-shell action of the D5-brane 
with the free energy of the theory. Let us introduce a cut-off at infinity, rmax, The wick rotated 
on-shell action is given by: 



S = Nf^AnVsR^H^/^ dffWl + l'\ 1 + ^ 

^ 9s J V ir^ + l'' 



(59) 



where V3 = J d?x and l{r) is the solution of equation (55 



expansion from equation (52), that the integral in equation 
large t^q^x- 



(f2+/2)2 

. It is easy to verify, using the 



59 ) has the following behavior at 



/ drrWl + l'\ 1 + ^-^ = lrl,^ + o(—] 

J V (f2 + /2)2 3 — \r^,^J 



(60) 



It is important that in these coordinates the divergent term is independent of the field /, it is 
therefore possible to regularize the on-shell action by subtracting the free energy of the / = 
embedding. The resulting regularized expression for the free energy is: 



where 



/d5 = / dr 




(f2 + /2)2 



(61) 



(62) 



A plot of /d5 vs. |m| is presented in figure |6| The states from the lowest positive branch in 
figure [6] have the lowest free energy and correspond to the stable phase of the theory. Therefore 
there is a spontaneous breaking of the global SO {3) symmetry and the theory at m = develops 
a negative condensate proportional to — Ccr ~ —0.59. Note that only the absolute value of m 
corresponds to the bare mass of the fundamental fields. The states with negative m correspond 
to D5-brane embeddings that intercept the / = line in the / vs. f plane and as seen from 
figure |6] are unstable. It is to be expected that the meson spectrum of the theory in such a 
phase would contain tachyons based on an analogy with the meson spectrum of the D3/D7 
system studied in ref. [T2|. Before we proceed with the analysis of the meson spectrum of the 
theory let us write an expression for the condensate of the theory {ipip) oc — Ccr = R^Hccr at 
zero bare quark mass. The coefficient of proportionality is given by [33j: 



iiT a — Cc 



lQTi^a'^^c,,R\H/2Txa'). 
9s 



(63) 
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Figure 6: The states corresponding to tlie lowest positive branch of the plot in figure [5] have the 
lowest free energy and thus correspond to the stable phase of the theory. 



Note that the condensate is proportional to the magnitude of the magnetic field Hjli^ol . In 
order to check the consistency of our numerical analysis and to calculate more accurately the 
constant Ccr we have calculated the value of Ccr for a range of H having set i? = 1. The resulting 
plot is presented in figure [Tj The solid (black) line corresponds to the linear fit Ccr ~ 0.586iJ 
therefore we have Cr-r ~ 0.586. 




Figure 7: A plot of q 
Ccr ~ 0.586i7. 



vs. 



H for R = 1. The solid (black) line corresponds to the linear fit 



3.3 Meson spectrum and pseudo-Golstone bosons 

In this section we will analyze the normal modes of the D5-brane. These describe fluctuations 
of the spinor bilinear in the dual gauge theory and hence their spectrum is the spectrum of 
the light meson-like excitations of the gauge theory. We focus our analysis on the normal 
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modes corresponding to the Goldstone bosons (which we label as pions here for simplicity) of 
the spontaneously broken 5*0(3) symmetry and study their spectrum as a function of the bare 
quark mass rriq. Our study shows that the external magnetic field splits the degeneracy of the 
meson spectrum and gives mass to one of the pions of the theory. It also modifies the standard 
oc m GMOR relation for the remaining Goldstone mode to a linear relation M^r oc m. We 
will show that these results are in accord with the behavior expected from the effective chiral 
Lagrangian of the theory. 

In order to study the light meson spectrum of the theory we look for the quadratic fluctuations 
of the D5-brane embedding along the transverse directions parametrized by /, ip, 4>- To this end 
we expand: 

1 = 1 + 27ia'6l; = 2T:a'5ip] = 27ra'50 , (64) 
in the action (50) and leave only terms of order {2Tia'Y. Note that fluctuations of the U{1) 



gauge field F^p of the D5-brane will also contribute to the expansion. There is also an additional 
contribution from the Wess-Zumino term of the D5-brane's action: 

5'wz = Nf^i^ /" ^ [Cp A e^] ; T = B + 2Txa'F . (65) 

For the anzatz that we are considering, the relevant term is: 

^wz = iV//i5 J BA P[C4 , (66) 

Me 

where P[C4] is the pull-back of the magnetic dual, C4, to the background C4 R-R form. For 
the particular parameterization of 5*5 considered here, it is given by: 

1 4r2/2 

C4 = — — —R^sinipildr - rdl) A dVL2 A dcj) . (67) 

Qs [r^ + vy 

After some long but straightforward calculations we get the following action for the quadratic 
fluctuations along I: 



5f , (6J 
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and along and ip: 

^f^M °^ \^/^S''^{G^^dJijdp5ij + G^^dJ(t>dp8<t>) , (69) 
oc {cosa)PH6ipdoS(f) . 

Here Eajj is the pull-back of the generalized metric on the classical D5-brane embedding: 

E^p = d^X^^dpX^iG,, + B^,) , (70) 

and we have defined S"^ and J"'^ as the symmetric and anti-symmetric elements of the inverse 
generalized metric E°'^: 

^ap ^ gaf3 ^ jaf3 _ ^7^^, 

The determinant E and the function K = P are given by: 



= {cos a)r^VTTJ^Jl + Ij;^—^ = g{r) cos a , (72) 

P = ^i^si-^' - ■ (73) 

As one can see, the fiuctuations along ip and decouple from the fiuctuations along / and the 
fluctuations of the gauge field A^. Since we are interested in the pseudo-Goldstone modes of 
the dual theory we will focus on the fiuctuations along ip and (p. The equations of motion 



derived from the quadratic action (69) are the following: 



where 



° = -a? + 7fi- (75) 

Note that the background magnetic field breaks the 5*0(1,2) Lorentz symmetry to 5*0(2), 
which manifests itself in the modified laplacian (75). Next we consider a plane- wave ansatz: 

S(j) = e^('^*-^^)r/i(r); Sip = e^(^*"^^)r72(r) , (76) 

now using the anzatz ( [76) ) we get: 
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The equations of motion in (77) can be decoupled by the definition ri± = rji ± irj2. The result 
is: 

g{r)f ,\ g{r)RH^ '^2 



(r2+P)2 



Because of the broken Lorentz symmetry, the 1 + 2 dimensional mass = — /c^ (depends 
on the choice of frame. We can define the spectrum of excitations as the rest energy (consider 
the frame with k = 0) and as we shall observe, the spectrum is discrete. Furthermore just as in 
the D3 /D7 case there is a Zeeman splitting of the spectrum due to the external magnetic field. 
Interestingly, at low energy the splitting is breaking the degeneracy of the lowest energy state 
and as a result there is only one pseudo-Goldstone boson. Note that this is not in contradiction 
with the Goldstone theorem because there is no Lorentz symmetry. This opens the possibility 
of having two types of Goldstone modes: type I and type II satisfying odd and even dispersion 
relations correspondingly. In this case there is a modified counting rule (ref. [2S], see also 
ref. [31]) which states that the number of GBs of type I plus twice the number of GBs of type 
II is greater than or equal to the number of broken generators. As we are going to show below 
the single Goldstone mode that we see satisfies a quadratic dispersion relation (hence is type 
II) and the modified counting rule is not violated. Note also that for large bare masses m (and 
correspondingly large values of /) the term proportional to the magnetic field is suppressed and 
the meson spectrum should approximate to the result for the pure AdSs x 5*^ space-time case 
studied in refs. [2T1 [22], where the authors obtained the following relation: 

o^n = |?V(n + l/2)(n + 3/2) , (79) 

between the eigenvalue of the n*^ excited state uJn and the bare mass m. 



In order to obtain the meson spectrum, we numerically solve the equations of motion (78) in 
the rest frame {k = 0). The quantization condition for the spectrum comes from imposing 
regularity at infinity. More precisely we require that ri± ~ 1/r at infinity (r oo). The results 
are summarized in figure [s] Just as in the D3/D7 case we have defined the dimensionless 
quantities m = m/R\/H and u = ujR/y/H. As one can see from figure ^ for large fh the 



spectrum asymptotes to that of pure AdSs x S"^, given by equation (79). The Zeeman splitting 
of the spectrum is also evident. It is interesting that as a result of the splitting of the ground 
state there is only a single pseudo-Goldstone mode. Furthermore, as can be seen from figure |9l 
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for small bare masses instead of the usual Gell-Mann-Oakes-Renner relation we obtain a linear 
dependence ~ m. As we will show in the next subsection the slope is given by the relation: 

4c 

= —m ^ 0.736m . (80) 




Figure 8: The meson spectrum of the first three excited states is plotted. There is Zeeman sphtting 
of the spectrum and the existence of a mass gap at m = as weh as a single Goldstone boson mode. 



For large fh the spectrum asymptotes to that of zero magnetic field given by equation ( 79 ) (straight 
lines). 

It is also interesting to study the dispersion relation of the Goldstone mode. Since we have 
broken Lorentz symmetry and observe only one pseudo- Goldstone mode (which is only half 
the number of broken generators) we anticipate a quadratic dispersion relation (see refs. [31] 
and [32] for discussion). 

In order to obtain the dispersion relation of the Goldstone mode we numerically solve equations 
(78) at very small bare mass fh 0.0007 and for a range of small momenta k = kR/\fH. The 



result is presented in figure 10 There is indeed a quadratic dispersion relation. As we are going 
to show, the dispersion relation is given by: 



r.2 



UJ 



'jk + 



TT 



CcrTn 



where: 



7 



dr- 



P/Vl + /'2 



+ P)Jl + (f^ + 1 



2^2 



For m ^ 0.0007 the relation (81) is given by: 



=.2 



UJ 



0.232A; +0.000515, 



(81) 



(82) 



^3) 
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Figure 9: Plot of the spectrum of the ground state from figure [s] for small bare masses. The dashed 
line corresponds to the linear behavior from equation ( 80 ) . 



and is represented by the fitted curve in figure 10 
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Figure 10: Plot of the dispersion relation of the pseudo-Goldstone mode for m ~ 0.0007. The 



parabolic fit corresponds to equation (83). 



In the next subsections we will obtain the effective 1 + 2 dimensional chiral action describing 
the pseudo-Goldstone mode and argue that in the limit c<j — > it is identical to the action de- 
scribing magnon excitations in a ferromagnet [11]. Furthermore we will show that the observed 



dispersion relation (81) is in agreement with the dispersion relation of magnons in an external 



magnetic field. Note that in order to make the analogy with a ferromagnet, one needs to iden- 
tify the bare mass with the external magnetic field acting on the ferromagnet. The reason is 
that these both correspond to the small parameter that explicitly breaks the global symmetry. 
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3.3.1 Low energy dispersion relation 



In order to obtain the dispersion relation for the pseudo-Goldstone mode we will analyze the 



first equation in (78) in the spirit of the analysis performed in section 2.2.1 for the D3/D7 
system. To begin with let us consider the limit of small uj thus leaving only the linear potential 
term m. Co. In view of the observed quadratic dispersion relation (81) we will also keep the k"^ 

g{r)RH^ 



term in equation ( 78 ) . 



dr 



l + /'2 



^4 



ujPH + 



V2 + Z2)2 + RiH2 



It is convenient to define the following variables: 

g{r)f 



l + /'2 



Then equation (84) can be written as: 

e 



^ — ^ 



02 



. 



(84) 



(85) 



(86) 



(r2 + /2)2 + ji^H^ 

Where the overdots represent derivatives with respect to r. Now if we take the limit m 
we have that uj —>■ and k —>■ and obtain that: 



uj=0 



e 



(87) 



is a solution to equation (86). Our next step is to consider small m and expand: 

^ = 9 + 5^; = 6 + 56, 

where the variations 6C, and 6Q are vanishing in the m — >■ limit. Then, to leading order in m 
(keeping in mind that a; ~ m and k'^ ~ m) we obtain: 

gir)RH' - 







LOPH 



(r2 + /2)2 + J14H2 











(89) 



Now we multiply equation (89) by and integrate along r. The result is: 



{eSi - 0(50 7 - (0^0 - 0(50) -ujH I drP{r) - ^R^VH-fk^ = 



(90) 



Using the definitions of 0, P{r) and ^ and requiring regularity at infinity for t]^, one can show 
that the first term in equation (|90| vanishes and that: 



(0(50 - 0(50) 



c6m 



drP{r) = -R^n/A 



(91) 
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and hence using the previous definitions, m = m/R^/H, c = c/R^H, uj = luR/^/H and 
k = kR/y/H, we obtain equation (81) which we duphcate below: 



UJ 



~2 4 

7fc H — Ccrm 



TT 



(92) 



In the next subsection we will derive the effective 1 + 2 dimensional action for the pseudo- 
Goldstone mode and show that to leading order it is in to one correspondence with the ef- 
fective action describing magnon excitations in a ferromagnet corresponding to the S0{3) 
5*0(2) spontaneous symmetry breaking by spontaneous magnetization [11]. We will relate the 
fermionic condensate c to the spontaneous magnetization of the ferromagnet and the bare mass 
to the external magnetic field and show that the dispersion relation (92) is in exact agreement 
with that of magnons. 



3.3.2 Effective chiral Lagrangian 

In order to obtain the 1 + 2 dimensional effective action describing the pseudo-Goldstone mode 



we consider the 1 + 5 dimensional action (69) for a classical embedding in the vicinity of the 
critical embedding, namely that embedding corresponding to a very small bare mass m. Now 
let us consider the following ansatze for the fields 50 and Sip: 

= e(;:y^i(^) ' = ' 

Since we are close to the critical embedding we will consider the same expansion as in equa- 



tion (88): 



^. = e + , i = 1 or 2 ; e = e + 59 . (94) 
By definition it follows that as m — * 0, 6C,i and 6Q vanish. Then to leading order we have that: 

dM=^[{mi-^S^i) + {^SQ-Q6Q)]xi{x) ■ d^6^ = d^xi{x) ; = 0, 1, 2 , (95) 

drSiP = ^[{m2-m2) + {^SQ-Q6&)]x2t{x) ; 9^5^ = 9^X2 (x) ; /x = 0,l,2. 



Now we integrate equation (69) along r from 0, oo and along the internal unit sphere fl2- The 
interesting term is the part of the kinetic term involving derivatives along r. After integration 
by parts it boils down to a mass term for the 1 + 2 dimensional fields xi? X2- Explicitly: 

drdn2^^^dr6(l)dM = -[ drdn^^d, (f^^r^^) H = (96) 



-47r[(e56 - 05^i) + (050 - 050)] 



^1 2 1 2 

2^^ 2^^ 
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Here we have used the same arguments as in equation (90). It is clear that one can perform 
an analogous calculation for the term involving drSif). The rest of the terms are dealt with 
straightforwardly by integrating along r. The resulting action is: 



I i -fdoX*doX - ^-fda*da - f^'^ixdoX* - X*doX) + '^{^^)oX*X \ , (97) 



where we have defined a complex scalar field x = Xi + ^X2- The constants in the effective action 
are given by: 

„ oo oo 

nil AT f,g{r)RH' . ^\ AT r ^ g{r)RH' 



m 



The effective action (97) is very similar to the one considered in ref. [3T], where the author 
studied Goldstone bosons in linear sigma models with chemical potential, only we have further 
broken the Lorentz symmetry by the introduction of an external magnetic field. The peculiar 
feature of this effective Lagrangian is the single time derivative term that is responsible for the 
unusual quadratic dispersion relation. In ref. [31] the authors have shown that the number of 
goldstone modes with quadratic dispersion relation is half the number of broken generators. 
This is exactly what we observe here (two broken generators but only a single goldstone mode). 

On the other hand, for the pseudo-Goldstone mode uj ^ and to leading order the effective 
action (97) can be written as: 

= j d^xlhle^^d^U'^U' - J^ph^U''^ - ^F^diU'^d.U''^ ; a = lor2 , (99) 



where 



t,» = (2W),"; E = 2,; = ^ ; ,3h ^ ^ c„,n'-^ . (100) 



Equation ( 99 ) corresponds to the effective action describing a ferromagnet in an external mag- 
netic field h with a spontaneous magnetization S and a magnetic coupling (3 [11]. Here 
U = {U^,U'^,U^) is a unit vector corresponding to the direction of the spontaneous mag- 
netization of the ferromagnet, and the action (99) is describing quadratic fluctuations of the 
magnetization near the ordered state = 1. The fact that the effective external magnetic field 
h is proportional to the bare mass m is to be expected since in both descriptions these are the 
small parameters that reduce the exact global symmetry to an approximate one by coupling to 
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the corresponding order parameters (magnetization an d qua rk condensate corresponding). If 
one takes into account the re-definitions from equation (100) and the definitions of u) and k, it 
is straightforward to check that the dispersion relation of ferromagnetic spin waves 

uj = -fk^ + Ph; 7=^, (101) 



is exactly that of equation (81). 

Of course many different microscopic systems exhibit the same low energy behavior and hence 
are described by the same effective Lagrangian. Furthermore the fact that the mass generation 
process in the D3/D5 system is associated to precisely the same global symmetry breaking 
pattern (50(3) —>■ 5*0(2)) as the transition from paramagnetic to ferromagnetic phase is also 
very suggestive. However it is the peculiar single derivative term coming from the Wess- 
Zumino contribution that is responsible for the observed dispersion relation. Here we have 
not investigated how far one can go in describing properties of ferromagnets using the D3/D5 
set up. It would be interesting to study interaction terms, which would require expanding the 
effective action beyond quadratic order. In any case it is somewhat satisfying that the dispersion 
relation of pseudo-Goldstone modes can be related to a real condensed matter phenomenon 
such as magnon spin waves. 



4 Universal properties of magnetic catalysis in Dp/Dq 
systems 

In this section we focus on some universal features of the mechanism of spontaneous symmetry 
breaking in an external magnetic field in the context of the Dp/Dq system. In particular we 
explore the observed discrete self-similar behavior of the theory in the vicinity of the trivial 
embedding / = corresponding to the non-symmetry breaking phase. As we mentioned in 
Section 2 and Section 3, for the D3/D7 and the D3/D5 systems, this embedding is unstable 
and the instability is manifested as a multi-valuedness of the equation of state in the condensate 
versus bare mass plane (m, c) seeded by a logarithmic spiral structure (see figure [2] and figure |5|. 
This spiral structure has been explored in details in ref. [5] in the case of the D3/D7 set up. It 
has been shown that the spectrum of meson-like excitations also exhibits discrete self-similar 
structure in the tachyonic sector of the theory. It is interesting that the same structure appears 
for other first order phase transition in the Dp/Dq set up, such as the meson melting |33j and 
electrically driven insulator/conductor phase transitions [M], [35]. In ref. [33j it was pointed 
out that the critical exponents (or more appropriately "scaling exponents" ) characterizing the 
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logarithmic structure exhibit universal properties and depend only on the dimension of the 
internal 5" sphere wrapped by the Dq-brane. A similar analysis was performed in ref. [35] for 
the case of electrically and R-charge chemical potential driven phase transitions and another set 
of "scaling exponents" was obtained. Here we will extend the analysis of the spiral structure in 
the case of the D3/D7 system performed in ref. [5] to the general case of Dp/Dq systems T-dual 
to the D3/D7 intersection and will show that the corresponding scaling exponents guarantee the 
existence of a discrete self-similar behavior in all Dp/Dq systems of potential phenomeno logical 
interest. This suggests the universal role of the external magnetic field as a strong catalyst of 
mass generation. 

To begin with, let us consider the zero temperature Dp-brane solution, given by: 

ds^ = KP i-dt^ + dxl \ + k} {du^ + u^d^l_p) , (102) 



i=l 



9: 



i-f^p 1 <-'01...p — J^n 



where Kp{u) = {R/uy~P and R is a. length scale (the AdS radius in the p = 3 case). Now 
if we introduce a Dg-brane probe having d common space-like directions with the Dp-brane, 
wrapping an internal S*" C 5*^"^ and extended along the holographic coordinate u, we will intro- 
duce fundamental matter to the dual gauge theory that propagates along a ((i+ l)-dimensional 
defect. 

Next we parameterize the transverse 9 — p plane du^ + u'^dfll_p by: 

dp^ + dL^ + p'^dnl + L^dn'^j_p_^ , (103) 

where dfl"^ is the metric on a unit radius m-sphere and p'^ + L'^ = u^. We also introduce an 
external magnetic field H /2'Ka\ corresponding to the Fp_i p component of the field strength 
tensor, by fixing a constant 5-field in the (xp_i,Xp) plane: 

B{2) = Hdxp_i A dxp . (104) 

Then the DBI part of the Lagrangian governing the classical embedding of the probe is given 

C oc e-yr^ = ^^p"yrT77^Jl+ ^ff"L • (105) 

^ 9s y (p^ + L^) 2 



^We consider only systems T-dual to the D3/D7 one, which imposes the constraint p — d+n+1 — A. 
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The equation of motion for the classical Dg-brane embedding is given by: 

mW-P \ 7-p p"Vl + L'2 LH^R^~P 



(p2 + L2)V 



(p2 + 



9-p 



(106) 



For large p ^ oo the second term in equation (106) vanishes and the solution L(p) has the 
asymptotic behavior: 



Up) = m H + 



P 



(107) 



which encodes [TS] the bare quark mass niq = m/27ra' and the quark bilinear condensate 
(ipip) oc — c of the dual gauge theory. 



It is also clear that the equation of motion (106) has a trivial solution L(p) = 0, which preserves 



the rotational symmetry in the 8 — p — n plane transverse to both the Dp and Dg-branes. This 
solution has zero bare quark mass and corresponds to the non-symmetry breaking phase of the 
dual gauge theory. The solutions in the vicinity of L = are unstable and correspond to the 
interior of the spiral structure that we are studying. In order to obtain the scaling exponents 
characterizing the spiral we will zoom in on the region close to the origin of the (p, L) plane. 
We first introduce dimensionless variables via: 



p = pRH''-p] L = LRH''~p] m = mRH''-p] c = cR"-H''-p; 



and now rescale: 



p = Xp ] L = XL . 



In the limit A — ^ equation (|106) becomes: 

L' 



P 



7 — p 



(p2 + U 



7-P 



1 + 



1 + L'2 



p"L 



(p2 + L2)^ 



. 



:io8) 



(109) 



(110) 



The solutions to equation (110) have the scaling property that if L{p) is a solution, then so is 



j^L{pp). In order to explore the vicinity of the critical L = solution we define L = + C(p) 
and linearize with respect to (, the result is: 



d.ip-'-^C) + ^P"-^C = . 



Ill) 



Next we look for solutions of equation (111) of the form C = ■ The quadratic equation for u 
that we obtain is: 

2z/2 + (n + d-6)i/+(n-c/ + 4) = . (112) 
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We have used the constraint p = 3 + d — n. Now in order to have a logarithmic spiral (which 
seeds the multi-valuedness of the equation state) we need to have two complex roots. The 
condition for that is: 



{n + d-6f <8{n-d + 4:) . 



;il3) 



Note that in order to be able to turn on a magnetic field we need d > 2. In addition we are 
not interested in theories with d > 3. It is then easy to check that for all possible values of n 



(clearly n < 5) the condition (113) is satisfied. Then the roots of equation (112) i'± are given 
by: 



n + d — 6 



> 



3 
4 



(n 



d + A)-{n + d~Qy . (114) 



The inequality in the second formula in equation (114) is saturated for the minimum possible 
values {n,d) = (1,2). The most general solution of equation (111) can then be written as: 

1 



C(P) 



P 



rn.d 



(ylcos(a„,dlnp) + i?sin(Q;„,rflnp)) 



(115) 



Now the scaling property of equation (110) suggests the following transformation of the param- 
eters {A,B) under re-scaling of the initial condition 1/(0) = Lq — > -Lq: 



A' 
B' 



cos (a„ In /i) sin (a^ In yu) 



;ii6) 



^r„+i sin (a„ In /i) cos(a„ln/i)^ 

For a fixed choice of the parameters A and i?, the parameters {A\B') describe a logarithmic 
spiral, whose step and periodicity are set by the real and imaginary parts of the critical/scaling 
exponents r„ ^ and an^- Note that from the inequality in equation (114) it follows that r„ ^+1 > 
I > and hence the spiral is revolving as one scales away from the critical L = solution. 

This self-similar structure of the embeddings near the critical solution L = in our zoomed in 
region parameterized by (p, L) is transferred by a linear transformation to the structure of the 
solutions in the (m, c) parameter space. The parameters corresponding to the critical L = 
embedding are given by (0, 0). Then sufficiently close to the critical embedding we can expand: 



M 



+ 0{A^,B^,AB) . 



;il7) 



The constant matrix M depends on the properties of the system. Generically it should be 
invertible (numerically we have verified that this is the case) and therefore in the vicinity of the 
parameter space close to the critical embedding (m, c) there is a discrete self-similar structure 
determined by the transformation: 

1 



/i 



r-n+l 



M 



COS {an In /i) sin (q;„ In n) 
— sin (a„ In p) cos(a„ In p) 



M- 



(118) 
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Note that the hnear map corresponding to the constant matrix M would rotate, stretch and/or 



shrink (along the different axes) the spiral defined via the transformation (116). However the 



overall shape of the curve defined via equation (118) still remains a spiral revolving around 



the origin of the (m, c) plane (see figure |5] for the case of the D3/D7 system). This suggests 
that the state corresponding to the center of the spiral (the L = solution) is unstable and 
hence there is a dynamical mass generation in the theory. (The stable state at zero bare quark 
mass has a non-zero condensate) Therefore we learn that for all Dp/Dq systems T-dual to the 
D3/D7 intersection (and with > 2 so that a magnetic field can be switched on) the effect of 
the magnetic field is to break a global internal symmetry and generate a dynamical mass. 

To conclude this discussion we will provide a numerical check of the consistency of our analysis. 
To this end we consider the separation of the Dq and Dp branes at Lj„ = -L(O) (note that 
is proportional to the dynamically generated quark mass). Now if we start from some and 
transform to Lj„ = we can solve for ^ and generate a parametric plot of fh / [LinY"''^^^ 



vs. an,d^og Lin /27T. The transformation (118) requires that the resulting plot be an harmonic 
function of unit period. For the particular case of the D3/D5 system we have r2,2 = ~l/2 and 
«2,2 = The corresponding plot is presented in figure 11 For sufficiently small Lj„ the 

plot is indeed an harmonic function of unit period. 




log(Lin) 



-1.5 



Figure 11: A plot of ^logL, 
function of unit period. 



vs. 



~ 1 /2 ~ 

rh/L-l^ . For sufficiently small Lin the curve is an harmonic 
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5 Conclusion 



In this paper we have investigated the properties of strongly coupled, large gauge theories 
in the presence of an external magnetic field. Both in three and four dimensions, the holo- 
graphic approach reproduces the behavior expected from classical field theory arguments and 
the magnetic catalysis of global symmetry breaking is shown to be a universal feature of a 
family of strongly coupled gauge theories. This is another success for the internal consistency 
of the gauge/gravity duality. 

As has now become commonplace in such studies, the spectrum of mesonic states is directly 
identifiable by numerical methods. Moreover, in the present study we are able to make a large 
number of statements analytically by studying the systems in the chiral limit. Indeed in this 
regime, in which the global symmetry is broken dynamically due to the presence of the magnetic 
field the results are almost completely tractable analytically. 

In 3+1 dimensions we are able to show that the Gell-Mann-Oakes-Renner relation holds exactly 
and in this setting where a subgroup of the special relativistic transformations remains unbroken 
in the presence of the magnetic field, a relativistic dispersion relation is indeed recovered. In 
the 2 + 1 dimensional case however, all trace of the boost invariance is lost once the magnetic 
field is turned on. The counting of Goldstone modes then becomes more subtle but we are 
able to show that this holographic setup gives the correct number of massless modes expected 
from the non-relativistic Goldstone counting rules. In addition we are able to show that these 
modes obey a quadratic dispersion relation, in contrast to the relativistic case in one spatial 
dimension higher. 

In the present investigation we have focused on the low energy chiral lagrangian, calculated up 
to quadratic order in the goldstone mode excitations. In the 2 + 1 dimensional setting we were 
able to show explicitly that at this order our system reproduces the low energy behavior of 
spin-wave excitations in a ferromagnet. The matching of symmetry breaking patterns between 
the two systems is the root of this equality. It seems unlikely however that such an agreement 
will hold to higher order. It would certainly be possible to investigate the systems discussed 
here at higher order in the low energy degrees of freedom and this would be an interesting 
direction for future work. 

The study of flavor degrees of freedom in a wide range of conditions is certainly well worth 
pursuing further using the methods illustrated in this paper. The magnitude of the phase- 
space involved with studying flavor, both abelian and non-abelian, in the presence of external 
electric fields and magnetic fields, temperature and finite chemical potential in a variety of 
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dimensions means that there is surely plenty more to be found in even the simplest systems. 

In the case of 2 + 1 dimensions, the understanding of gauge theories in the presence of low 
temperature and high magnetic field is of much interest. In particular this is one area which 
may be accessible experimentally and one may be able to gain insight into such intriguing 
phenomena as the quantum hall effect. It seems very likely that we will make further inroads 
into understanding such effects using the AdS/CFT correspondence in the not too distant 
future. 

The diversity of phenomena that we can investigate using holographic techniques is clearly 
far larger than was expected in the early days of the AdS/CFT conjecture. The prospects of 
obtaining deep insight into such fascinating systems as non-conventional superconductors, the 
quantum hall effect and strongly coupled plasmas are real and exciting and the community 
continues to make progress in these directions. 
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Appendices 



A Low energy effective action for <l> 



Let us begin with the Lagrangian for the quadratic fluctuations equation (11) 



9s ^ 



(119) 



9s 



Caa = -{27ra'y^y^\\gS'"''S'''FatFa'b' • 
9s 4 

In order to obtain an effective action for the fluctuations along $ we will integrate out the 
fluctuations of the gauge field and in particular the Aq, Ai components. In more detail the 



contribution from the last two equations in equation (119) can be written as: 

1 



.11 ^aa'( 

+HdpK^Foi 



Caa + Ua cc -gS''S^''■i-^aAo^a'Ao + ^aA^^a^A^)--g{S'')\^oAo-^,A,f (120) 



To integrate out the Aq, Ai components of the gauge field we simply obtain the equations of 
motion for Aq and Ai and substitute them into the action. The equations of motion are: 



daigS^'S'^'^'daAQ) + giS^'fOoidQAQ - ^i^i) + Hd.Kdi^ = 
daigS^'S'^'^'d^A,) + giS'^fd^idQAQ - ^^A,) + HdpKdo^ = 



(121) 



which can be written as: 



da (^5"5""'9,Foi) - ifa,irn(i,i)<i> = ; □(!,!) = -dl + dl , 
da {^gS^'S^'^'daidQAQ - 9iAi)) - (7(5'')'n(i,i)(9oAo - d,A,) = . 



(122) 



Substituting back into the action ( 120 ) and integrating by parts leads to: 



Caa + cc --M9p($Foi) . 



The equation of motion for Fqi can be written as: 

^P^oi = ™°(i.i)<^' - ^ / dpgiS'r^FQ, , 

1 1 



(123) 



(124) 
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where is defined in equation (22). Now we substitute into the action (123) and denote by 
$(a;) the dimensionally reduced field $, to obtain: 



9s 



2 *f 



-$n(i_i)<l> + . . . 



(125) 



where we have ignored higher order derivatives terms. Combining this with the dimensionally 
reduced term we obtain the result: 



(126) 



Where u and ^ are defined in equation (22). In order to obtain a mass term for the dimen- 
sionally reduced field $ we have to take into account the radial dependence of the field Our 
analysis from Section 2.2.1 suggests that we should consider the following ansatz: 



<|.(p,x) = ^|.(x) 



(127) 



where ^ is defined in equation (22) and we require that for the spontaneous symmetry breaking 
classical embedding (denoted by Zq) we have that ipllo = ^Ilq = ^- Then if we consider 
embeddings in the vicinity of Lq corresponding to small bare quark mass 6m we can expand: 



iP = \iJ + Stp ; ^{p,x) 



1 + 6 



(128) 



Now if we demand that as 6m ^ we have small momenta and a small mass term (which 
vanish at the critical embedding) to leading order we still have the expression from equation 
(126) plus some small mass term involving derivatives along p: 



2 ^ 



$2 + 



(129) 



where we have integrated by parts the last term and the dots represent higher derivatives terms 
and other sub-leading terms. Now it is straightforward to integrate along the unit and the 
radial coordinate p. Let us provide some more details in the integration of the mass term: 



dpdp 



-2c6m^ 



(130) 



Then for the final form of the effective action one obtains: 



S, 



eS 



U d'xi |-(9o*)' + + 7l(S2t)' + (93*) 



^m,4n + .... (131) 
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where: 



Af ={27ra' fNf^TT^ / dp 







2 AAf Sm 



/' = 7T^^-7^^; "^. = 7^. (132) 



^2 J ' (27ra')^' 27ra' 



.^/*(5tV/^p(^t^ + ^); (W^-pSp^. (133) 







One can see that this is the most general quadratic action consistent with the S0{1, 1) x S0{2) 
space-time symmetry. Furthermore the exphcit form of the mass term is in accord with the 



Gell-Mann-Oakes-Renner relation (37). To obtain the expression for provided in equation 
(132) one needs to consider the strict rrig — ^ limit and use that in this limit $ = 0/(27ra'). 
Next since corresponds to rotations in the transverse plane and is thus the angle of chiral 



rotation [18], the normalization of the kinetic term in the effective action (131) is given by 
A/" = (27r«')2/^/4. The last relation determines in terms of J\f. 
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